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rn \ ABSTRACT 

■ We have discovered T Tauri stars which show starthng spectral variabil- 

ity between observations seperated by 20 years. In spectra pubhshed by 
iBouvier fc Appenzeller fl992t) these objects showed very weak Ha emission, broad 
Call absorption and so called "composite spectra", where the spectral type inferred 
from the blue region is earlier than that inferred from the red. We present here new 
spectroscopy which shows that all four stars now exhibit strong Ha emission, narrow 
CaH emission and a spectral type which is consistent at all wavelengths. 

Dh' We propose a scheme to understand these changes whereby the composite spectra 

Q . of these stars can be explained by a period of active accretion onto the central, young 

' star. In this scheme the composite spectrum consists of a contribution from the stellar 

^ \ photosphere and a contribution from a hot, optically thick, accretion component. The 

■ optically thick nature of the accretion flow explains the weakness of the Ha emission 

during this phase. Within this scheme, the change to a single spectral type at all 
wavelengths and emergence of strong Ha emission are consistent with the accretion 
columns becoming optically thin, as the accretion rate drops. There is a strong analogy 

i-j \ here with the dwarf novae class of interacting binaries, which show similar behaviour 

during the decline from outbursts of high mass-transfer rate. 

The most important consequence of this interpretation is that these objects bring 
into question the association of Weak-Line T Tauri stars (WTTs) with non-accreting 
or discless objects. In light of this result we consider the justification for this paradigm. 

Key words: accretion, accretion disks, stars:pre-main-sequence, instabilities, 
stars:dwarf novae, planetary systems: protoplanetary discs 



1 INTRODUCTION 

T Tauri stars are divided into classical T Tauri stars 
(CTTs) and weak T Tauri stars (WTTs) largely on the 
basis of their Hq emission. The cutoff is usually taken 
at lOA, with values abov e this representing a CTTs 
l|Appenzeller fc Mundt 19891) . It has long been argued that 
the difference between CTTs and WTTs is that the for- 
mer have discs, and the latter do not. Indeed, there is some 
support for this paradigm, with Ha equivalent width corre- 
lating with other dis c indicators, such as B-V excess , or in- 
frared colour excess jHerbig 19981: ICabrit et al. 1990D . How- 
ever, many problems exist with such a scheme. The cutoff 
of lOA is somewhat arbitrary, and is less sensitive at earlier 



spectral types (iMartm 19971) . Also, recent studies have dis- 
covered WTTs w hich show infared excesses indica tive of cir- 
cumstellar discs llGregorio-Hetem fc Hetem 2002ft . Here we 
present evidence that at least some WTTs not only possess 
accretion discs, but are accreting from them at a rate which 
is comparable to the most rapidly accreting CTTs. 

An optical spectroscopic and photometric survey of a 
set of 47 stars associ ated with Einstein X-ray sour ces in the 
p Ophiuchus cloud (iBouvier fc Appenzeller 1992L hereafter 
BA92), identified 30 pre- main sequence stars. Four of these 
stars (ROXs 21, ROXs 47a, ROXs 2 and ROXs 3) were iden- 
tified as exhibiting "composite" spectra - the spectral type 
inferred from a blue spectral region near H/3 was systemat- 
ically earlier than that inferred from a red spectral region 
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surrounding the TiO band between 7050 and 7200 A. The 
spectral types inferred for these stars are listed in Tabled 
We note here that the composite spectrum effect was reliably 
detected in the case of ROXs 21, 47a and 3. The detection of 
the effect in ROXs 2 was less secure; although the blue spec- 
tral type of K3 was reliable, the lower resolution spectrum 
of this object rendered the red spectral type of MO uncer- 
tain. Our spectrum published here confirms the red spectral 
type of BA92, and hence we consider this star as a bona-fide 
composite spectrum object. BA92 also found a discrepancy 
in red and blue spectral types for ROXs 29, but the differ- 
ence in spectral type (K4 vs K6) was barely significant, and 
we do not consider this star a bona-fide composite spectrum 
object. 

BA92 put forward two suggestions for the origin of the 
composite spectra. 

• The systems are close binaries, with components of dif- 
ferent temperatures. 

• The systems consist of cool stars surrounded by accre- 
tion discs. The blue spectral region is dominated by luminos- 
ity either from the accretion itself, or from hot spots arising 
where the accretion stream hits the star. 

However, they were unable to determine which, if either, 
of these explana tions was the c ase. The binary hypothesis 
was examined by|Koreskoil99l, who performed speckle in- 
terferometry for ROXs 21 (SR 12), ROXs 47a (Do-Ar 51) 
and ROXs 3. He found that ROXs 21 and ROXs 47a were 
close binary stars with brightness ratios at 1.65 /im of 0.9 
and 0.33 respectively. However, he also found that co-eval 
binaries could not be responsible for the composite spec- 
tra effect. The problem arises as the two components must 
have comparable luminosities at 5550 A. The cool compo- 
nent must then be larger than the hot component, implying 
that the hot component is substantially older than the cool 
component. 

Hence, the nature of the objects discovered in BA92 
remains a mystery. In order to resolve this mystery, we 
revisited the objects, more than twenty years later, with 
medium-resolution spectroscopy in the blue and red regions. 
The spectra of all four objects show significant changes since 
the observations of BA92. All objects now show strong Ha 
emission and a diminished composite spectrum effect. The 
presence of such variability strongly argues against a binary 
origin for the phenomenon. Here we suggest that the com- 
posite spectrum phenomenon, and the variability witnessed, 
is caused by changes in the accretion state of the young 
stars. 

Section |5| describes the observations taken, whilst the 
results are shown in section |3 In section |1] we outline the 
scheme we propose to explain these results, and discuss its 
consequences and in section|K|we summarise our conclusions. 



2 OBSERVATIONS 

On the nights of 2002 Jan 27 and 28 we obtained 5960- 
7640 A spectra of ROXs 21, ROXs 47a, ROXs 2 and ROXs 3 
using the RGO spectrograph on the Anglo- Australian Tele- 
scope (AAT) at Siding Spring, Austrafia. The 1200R grat- 
ing, blaze-to-camera, in conjunction with the MITTL chip 
gave a resolution of 1.5A. The exposure time was 120 sees 
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Figure 1. 2dF blue spectra of the BA92 objects, obtained in 
August 2002. The spectra have been normalised and an offset has 
been applied to the y-axis of the spectra for clarity. 

for each object. The spectra were ext racted using the opti- 
mal extraction algorithms in FIGARO iCohen 19881) . Obser- 
vations of the A- type star SAO 184424 were used to correct 
the spectra for telluric features. Regular observations of arc 
and flatfield frames were used to calibrate the data. 

In addition, on the nights of 2002 August 30 and 2002 
September 1, we obtained 3770-4890A spectra of a sample 
of 23 of the 30 pre- main sequence object s oberved by BA9 2 
using the 2dF spectrograph on the AAT Jl ewis et al. 200211 . 
Included in the sample were the four composite spectra ob- 
jects. Each spectrum consists of 5x300 sec exposures. The 
data were taken in service mode by E. Corbett. The 1200B 
grating was used, giving a resolution of 2.2A. Arc and flat- 
field frames were obtained after each target exposure, along 
with offset sky frames (which may be used to calibrate fibre 
throughput). The data were reduced usi ng the 2dfdr dat a 
reduction package, with default settings iLewis et al. 200211 . 



3 RESULTS 

Figure0shows the 3770-4890A spectra of ROXs 21, ROXs 
47a, ROXs 2 and ROXs 3. Shown for comparison in FigureH 
are the original data from BA92. The variability is striking. 
First, whilst the later spectra of all four objects are charac- 
terised by emission lines of the Balmer series and Call, the 
original data show no Balmer emission, and strong, broad 
Call absorption. Only one object (ROXs 21) shows any sign 
of Call or Balmer emission in the ealier data. In all cases, 
the increase in emission line strength is significant and strik- 
ing, even when our data are degraded to the resolution and 
signal to noise of the BA92 spectra. Second, and most in- 
triguing, the prominent G-band seen in the original data 
at around 4300A is much weaker in the current data. The 
G-band is visible in dwarfs from mid-F to late-K type, and 
peaks in strength around G5V. 

With the weakening of the G-band, it is clear that the 
spectral type inferred from the blue spectral region is later 
in our data than in the data of BA92. In fact, BA92, used 
the continuum shape and G-band strength to infer spectral 
types of K2V-K4V for these objects. Comparison of our 
spectra and the K7/M0V spectrum in Figure |5| shows that 
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Figure 2. The original spectra of the BA92 objects, obtained 
between 1983 and 1986. The spectra have been normalised, and 
an offset has been applied to the y-axis of the spectra for clarity. 
Also shown is the spectrum of a K4 dwarf (bottom) and a K7/M0 
dwarf (top). 
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Figure 3. RGO red spectra of the BA92 objects, obtained in 
January 2002. The spectra have been normalised and an offset 
has been applied to the y-axis of the spectra for clarity. 



the spectral type inferred from the blue spectral region is 
now closer to late-K or early M. The implications of this are 
discussed in section |3] 

The red spectral region has also seen significant changes 
since the original data were taken. BA92 classified all four 
objects as weak T-Tauri stars on the basis of their Ha emis- 
sion. Figure O shows our recent AAT spectra. A glance at 
this figure shows that the Ha emission properties of the 
objects has changed markedly: the red spectra of all four 
objects are now dominated by strong Ha emission. Tabled 
shows the changes in Ha equivalent width between the two 
observations. In all cases the equivalent width of Ha emis- 
sion has dramatically increased. The change is large enough 
such that all four objects would now be classified as classi- 
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Table 1. Ha equivalent width variations between this paper and 
BA92. Also shown is the Ha full-width at 10% of maximum in- 
tensity for our data, and the spectral types inferred by BA92. The 
first entry is the spectral type inferred from the region bluewards 
of H/3, the second entry is the spectral type inferred from the 
region around Ha. 



Object 


BA92 


This paper 


This paper 


BA92 




EW (A) 


EW (A) 


FW10% (kms-i) 


Sp Type 


ROXs 21 


4.5 


17.6±0.2 


370±30 


K4/M2.5V 


ROXs 47a 


9.2 


24.4±0.2 


460±30 


K2/K7-M0V 


ROXs 2 


2.8 


14.3±0.2 


550±30 


K3/M0V 


ROXs 3 


2.3 


12.2±0.2 


500±30 


K3/M0V 



cal T Tauri stars, on the basis of their Ha emission. The Ha 
emission is also very broad. White & Basri (2003) show that 
the Ha line widths at 10% of peak intensity are often a bet- 
ter indicator of ongoing accretion than Ha equivalent width. 
They suggest that stars with 10% widths of > 270km s~^ are 
accreting, regardless of spectral type. All four objects show 
10% widths significantly in excess of this value (see TableQi. 

Despite the change in Ha emission properties, the spec- 
tral type inferred from the red spectral region remains un- 
changed with respect to that found by BA92. This removes 
any uncertainty that the objects we observed are the same 
objects observed by BA92. We now discuss the implications 
of these results. 



4 DISCUSSION 

4.1 The nature of composite spectra 

Any theory which seeks to explain the composite spectra 
observed in these stars must explain these key facts. 

• A composite spectrum always shows a hot spectral type 
in the blue, and a cooler spectral type in the red. 

• The phenomenon is transitory, and decays on a 
timescales of less than or of the order 20 years. 

• Coupled with the disappearance of the composite spec- 
trum, is the emergence of an emission line spectrum from 
CaH and the Balmer lines in the blue region, and a dra- 
matic rise in the equivalent width of Ha emission in the red 
region. 

Considering these facts, we can reject binarity as a cause 
for composite spectra; it is unlikely that the spectral type 
of a star can change by several subclasses in £ 20 years. A 
second possibility is that the composite spectra are caused 
by cool starspots on the surface of a hot young star. This, 
too, seems unlikely. On a spot-covered star, the overall light 
is normally dominated by the hot, bright, spot-free pho- 
tosphere. BA92 calculated that to produce the composite 
spectrum effect observed required 93% spot coverage. Given 
that the spectrum in the blue region has also now cooled, 
starspots would now have to cover the entire star to explain 
the observed changes. 
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The most promising scheme is that the composite spec- 
tra are caused by accretion. In this scheme the source for 
the hot, blue component is the accretion luminosity itself. 
Unlike the veiling luminosity typically seen in CTTs, which 
is featureless and fills in stellar absorption lines, we are sug- 
gesting that the accretion luminosity, arising from optically 
thick accretion columns or in hot spots on the stellar sur- 
face, itself exhibits hot spectral features, like the G-band at 
43OOA. The stellar spectrum is then completely hidden by 
this accretion luminosity. An optically thick accretion flow 
would explain the absence of strong Ha emission in this 
phase. An attractive property of the accretion model is that 
the observed spectral changes are easily understood as a fall 
in the accretion rate. First we consider the blue spectral re- 
gion. As the accretion rate drops, the accretion flow becomes 
optically thin, and the luminosity of the flow reduces dras- 
tically. The hot spectral features will naturally disappear 
at this point. The drop in accretion luminosity explains the 
change to a red continuum, and emergence of an emission 
line spectrum as the veiling caused by the luminosity disap- 
pears. In the red spectral region, the observed changes are 
simpler, because this region is relatively unaffected by the 
accretion emission. In this case, the main observed change 
is a rise in the strength of Ha emission as the accretion 
columns become optically thin. 

This behaviour is not without precedent. Dwarf novae 
are binary stars consisting of a low-mass red star transfer- 
ring mass, via an accretion disc to a white dwarf companion. 
They exhibit frequent outbursts, which are linked to an in- 
stability in the accretion disc which causes sudden increases 
in accretion rate. These objects are characterised in their 
quiescent state by a strong emission line spectrum. During 
outburst, however, high accretion rates mean the disc at- 
mosphere becomes optically thick, and an absorption line 
spect rum, arising from the disc ph otosphere, is observed 
fe. g. iMansperger fc Kaitchuck 1 990V The similarities be- 
tween this behaviour and the behaviour seen in the BA92 
objects are striking. 

The most obvious implication of this explanation is that 
here are four stars which may exhibit high accretion rates in 
the WTT phase. If the blue luminosity in the BA92 objects is 
dominated by accretion power, we can calculate the required 
accretion rates from the observed B and V magnitudes, after 
correction for reddening. We assume the accretion emission 
is a black body, with its temperature of 5200K taken to 
correspond to the blue spectral type observed by BA92. To 
reproduce the composite spectrum effect we scale the lumi- 
nosity of the accretion source until it contributes roughly 
twice the stellar luminosity in the V-band. We use the pho- 
tometry and reddening of BA92, although we note that the 
BA92 derived the reddening from the V-R colours of the 
objects. Since the intrinsic V-R would have been affected by 
accretion luminosity these reddening values are likely to be 
underestimated, implying that the accretion rates derived 
here are a lower limit. In order to match the observed B and 
V magnitudes implies the accretion source has a luminosity 
of 0.1 Lq. We assume a stellar mass of 0.5 Mq and an inner 
accretion disc radius of 5 R^, which corresponds to 3 stellar 
radii at an age of 1 Myr llD'Antona fc Mazzitelh 19971) . With 
these values a luminosity of 0.1 Lq imply accretion rates 
of at least lO~^M0yr~^, highe r than the avera ge value of 
lO"*M0yr-i found for CTTs (iHartmann 1998h . Ahhough 



our value of lO~^M0yr~^ is only a rough estimate, it is cer- 
tainly true that the composite spectrum objects were accret- 
ing at a rate comparable to that seen in CTTs. This would 
be a startling reversal of the paradigm in which strong Ha 
emission is taken as evidence of strong accretion onto the 
young star, and weak Ha emission denotes a non-accreting 
source. 



4.2 Near-infrared excesses 

It is pertinent to consider if there is other evidence for 
the existence of accretion discs in these stars. Circum- 
stellar discs can produce strong infrared emission signifl- 
cantly in excess of the stellar luminosity. Such excess in- 
frared emission is easily detected in a near-infrared colour- 
colour diagram (see lHaisch et al. 200ll . for example) . We col- 
lated JHK photometry for t he co mposite spectrum stars 
from a number o f sources l^ ouvicr & ADDonzcUcr 1993; 
'Wa her et al. 19 94: Rvdgrcn et al. 19761 : ICutri et al. 2003: 
IJensen fc Mathieu 1997 ) . The data were transformed to the 
UKIRT photom etric system using the transformations in 
ICarpenter (20011) and an average of each magnitude was 
taken. None of the four stars showed any evidence for a 
near infrared excess when placed on a JHK colour-colour 
diagram. L band photometry was also available for ROXs 
47a and 21. JHKL colour-colour diagrams are a much more 
sensitive disc indicator than JHK diagrams, often reveal- 
ing evidence for circumstellar discs in objects showing no 
JHK excess (see iHaisch et al. 200ll . for example). In this 
case, ROXs 47a showed a K-L excess of a few tenths of a 
magnitude, whilst ROXs 21 showed no K-L excess. The K-L 
excess i n ROXs 47a may, h owever, be due to a low-mass com- 
panion iSimon et al. 1987^ . Furthermore, we caution about 
over-interpretation of these data; T Tauri stars are intrinsi- 
cally variable in the near i nfrared, with varia bility up to 1 
magnitude not uncommon jEiroa et al. 2002^ . As the pho- 
tometry obtained above was not obtained simultaneously it 
seems plausible that variability could mask any near infrared 
excess present. 

Although there is a strong correlation between near- 
infared excess and other disc indicators, the absence of a 
near infrared excess does not nec essarily imply the a bsence 
of an accretion disc. For example. iRebull et al. f2002r) found 
several stars with very strong Ha emission, presumably aris- 
ing in accretion columns, which exhibited no near infrared 
excesses. This presumably means it is possible to maintain 
signiflcant accretion rates onto a young star through a disc 
which is optically thin in its inner parts. 

4.3 Time-dependant behaviour of the composite 
spectrum state. 

To better understand the physics of what might cause 
changes in accretion rate it is pertinent to consider the time- 
dependance of these changes. For example, how long does 
the high accretion state last? Is the transition between high 
and low states short? 

The BA92 sample consists of 30 X-ray selected pre-main 
sequence stars which also exceed a given optical flux limit. In 
their sample, 4 out of 30 objects show a composite spectrum 
effect. Naively, this suggests that ~10% of T Tauri stars are 
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in such a state (but see below). Twenty years later, our 2dF 
spectroscopy shows that all four objects are now in a low 
accretion rate state - this argues that the high state is likely 
to last less than 20 years, with long-duration high states of 
50 years or more being very unlikely. Furthermore, in the 
study of BA92, the composite spectrum effect was either 
absent or showed similar amounts of spectral type mismatch 
in all objects. In other words, no objects were observed in 
a transition between high and low states. This then implies 
that the transition time between the two states is short - 
statistical arguments suggest the transition time is unlikely 
to be longer than 0.5 years. 

The timescale of 0.5 years is interesting because it sug- 
gests a physical mechanism for the change in accretion state. 
In dwarf novae, the cause of the accretion rate variations is a 
thermal- viscous instability (see lLasota (200lh for a review). 
A similar thermal-viscous instability is believed to cause the 
FU Orionis phenomena in young stars. In dwarf novae, the 
accretion rate changes on the thermal timescale for the disc. 
The thermal timescale for the inner region of T Tauri discs 
is roughly half a year (at a radius of 0.1 AU). This suggests 
a thermal-viscous instability in the inner disc may be re- 
sponsible for the change in accretion rates between the two 
states. 

Interestingly, none of the nineteen other objects we ob- 
served from the BA92 sample underwent a change from a 
low to a high state - the spectral type inferred from our 
2dF spectra agreed with the classification of BA92, with the 
exception of the four composite spectrum objects^ . In other 
words, the sample of BA92 contained four high state objects 
and no low state objects. Although we have observed four 
high-low transitions, no low-high transitions have occurred 
within the BA92 sample. This suggests that the sample of 
BA92 was strongly biassed towards objects which were in 
a high state at the time of their observations, and that in 
fact the fraction of objects in a high state at any given time 
may be much lower than 10%. Only simultaneous blue/red 
spectroscopy of an unbiassed sample of pre-main sequence 
stars can resolve this question. 

4.4 Hq variability in T Tauri stars 

We note here that this is not the first evidence for large 
cha nges in Hq emission st rength in T Tauri stars. For exam- 
ple, Ijoergens et aL (200111 noted that the equivalent width 
of Ha emissi on in the T T auri star RX J1608. 6-3922 varied 
from 7-14A. iHerbig f 19981) noted that several of the young 
stars in the cluster IC 348 showed variable Ha emission, 
with some objects changing from 2A to > lOA in e quiva- 
lent width. In A-Ori, the survey of iDuerr et al. fl98^ cata- 
logued 82 stars with Ha emission, which they categorised as 
"s trong", "medium" and " weak". The region was revisited 
bv lDolan fc Mathieu ('2o'oih . who found that the Ha equiva- 
lent widths of the "strong" sample included lines as weak as 
GA, a fact they attribute to variability. Then there is the re- 
markable case of KH 15D - a weak T Tauri star which shows 

^ ROXs 4 did show a discrepancy with the BA92 spectral type of 
K4 being incompatible with our classification as early-M, however 
close inspection of the spectrum published in BA92 leads us to 
conclude that this star was mis-classified by the authors. 



eclipses thought to be caused by inner disc material. During 
these eclipses the Ha equivalent width ri ses from around 2A 
to nearly 30-50A jHamilton et al. 2002^ . 

Since none of these studies have incorporated spec- 
troscopy in the blue spectral region, it is not possible to 
say if these changes in Ha equivalent width are accompa- 
nied by the changes we show here. It is clear, however, that 
there are many examples in the literature where objects have 
changed their nature from classic to weak T Tauri stars and 
vice- versa. 



4.5 The link between Ha and discs in T Tauri 
stars 

It seems certain that strong (EW > 10 A) Ha emission 
is indicative of the presence of active accretion and hence, 
an accretion disc. Chromosphe ric emission from Ha rarely 
rises above lOA iMartin 1997 *1 . so an additional source of 
Ha emission must be found. Furthermore, in stars with Ha 
emission of EW > lOA, the stre ngth of Hq e mission is pro- 
portional to the B-V excess (e.g. IHerbig 199^ . T Tauri stars 
with strong Hq also show the infr ared excesses thoug ht to 
arise from the inner accretion disc iCabrit et al. 1990fl . 

However, it seems that the conclusion that WTTs lack 
discs is not correct. In fact, it seems certain that at least 
some WTTs do possess discs. For example, in the case of 
KH 15D, which is a WTTs, it is material in a circum- 
stellar disc that is believed to be r esponsible for ecl ipses 
in this system llHerbst et al. 2002ll . iBarv et al. ^2002^ de- 
tect molecular hydrogen, believed to originate in an ac- 
cretion disc, again around a WTTs. Also, the variability 
of Hq emission on timescales as short as twenty years ar- 
gues against its usefuleness as a disc indicator (it is un- 
likely that the discs themselves could be coming and going 
on these timescales). Furthermor e, in a sample of 27 WTTs, 
iGregorio-Hetem fc Hetem f20 02^ find that 60% of the ob- 
jects show infrared excesses which indicate a circumstellar 
disc. 

In at least some of these cases it is tempting to explain 
the failure of the accepted paradigm by stating that WTTs 
may or may not have discs, but that these discs are currently 
not accreting onto the central star. If our interpretation of 
the composite spectrum phenomenon is correct, however, we 
are forced to accept that not only are some WTTs accreting, 
but they are accreting at rates higher than, or comparable 
to most CTTs. 

The presence of active accretion discs amongst the 
WTTs population is an important result. It suggests that 
there is no fundamental difi'erence between CTTs and some 
WTTs, other than the current accretion rate. This might 
explain, amongst other things, why an evolutionary link 
between CTTs and WTTs has been difficult to establish. 
Furthermore, if the time-averaged accretion rate in T Tauri 
stars is higher than previously suspected, this might explain 
age spreads seen in some star forming regions. An accre- 
tion rate of lO~^M0yr~^, as seen in the BA92 objects is 
high enough to upset the thermal balance of a T Tauri 
star, moving it ba ck towards the birthline in a HR diagram 
jTout et al. 19991) . 
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5 CONCLUSIONS 

We have obtained new spec tra of the objects found to ex - 
hibit composite spectra by iBouvier fc Appenzeller f 19921) . 
The composite spectrum effect has diminished in the inter- 
vening years; the blue spectral type now agrees more closely 
with that inferred from the red spectral region. In addition, 
a series of emission lines has appeared, and the equivalent 
width of Hq emission has risen drastically. These spectral 
changes can be understood if the composite spectrum effect 
is a result of the blue spectrum being dominated by hot, opti- 
cally thick emission from the accretion flow. The luminosity 
in the blue then implies accretion rates of lO~'^M0yr~^. An 
important result of this explanation is that the paradigm re- 
lating Hof emission to accretion discs is incomplete; at least 
some WTTs both posess accretion discs, and are accreting 
from them at a high rate. 
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